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Abstract. We consider coherent single-electron dynamics in the one-dimensional 
nanostructure under resonant electromagnetic pulse. The structure is composed of two 
deep quantum wells positioned at the edges of structure and separated by a sequence 
of shallow internal wells. We show that complete electron transfer between the states 
localized in the edge wells through one of excited delocalized states can take place 
at discrete set of times provided that the pulse frequency matches one of resonant 
transition frequencies. The transfer time varies from several tens to several hundreds of 
picoseconds and depends on the structure and pulse parameters. The results obtained 
in this paper can be applied to the developments of the quantum networks used in 
quantum communications and/or quantum information processing. 
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1. Introduction 

The implementation of reliable connection between remote parts of quantum register 
remains one of the serious challenges on the way of realization of scalable quantum 
computations. One possible solution of this problem is to encode quantum information 
into the spin or charge state of single probe electron that can next be transferred from one 
qubit to another qubit with the help of an auxiliary quantum circuit. In the theoretical 
schemes of solid-state quantum information processing, such quantum networks based 
on regular arrays of specific quantum nodes (semiconductor quantum dots fT] , quantum 
wells [2j, quantum wires pi], donor atoms [^), exploiting voltage-controlled tunneling for 
coherent electron transport, were recently proposed. The transfer protocols developed 
in the most of those works employ the stimulated adiabatic passage techniques adopted 
from the quantum optics. These methods were shown to be quite promising due to their 
robustness against unwanted transitions from the transport state. The electron transfer 
time depends on the distance between the sender and target nodes as well as on the 
voltage pulse design. In Ref. [3] it was reported on the possibility of complete electron 
transfer between two electrostatically defined quantum dots connected through the 600 
nm long quantum wire. The calculated transfer time r ~ 50 ps is short enough to allow 
coherent manipulations with the nonlocal charge electronic state. However, since the 
voltage driving schemes are very sensitive to imperfections in the time delays and in 
the shapes of voltage pulses, the requirements on the pulse engineering remain severe. 
This makes the realization of such transfer protocols to be difficult experimental task. 
Other strategies of selective electron transfer utilize the strong and rapid electromagnetic 
modulation of the potential relief of nanostructures [5[ [6]. They rely upon the fact 
that the tunneling rate between neighboring nodes in the structure can be efficiently 
controlled by the field strength variation. At the same time, the strong fields acting on 
the auxiliary network will disturb the qubits in uncontrollable manner that may create 
insurmountable obstacle on the route of scalable quantum computer's building. 

In this paper, we propose an alternative mechanism for the selective electron 
transfer between two edge identical quantum wells of the quasi-one- dimensional 
multiwell nanostructure, where internal wells form the body of the transport circuit. 
The probe electron, occupying initially the ground state of one of the edge wells, can 
be brought into interaction with the nearest qubit. Thus, the probe electron becomes 
entangled with this qubit, and its spin or charge state is conditioned by the qubit state. 
Being transported between two edge wells, the electron carries the information about the 
state of a given qubit onto the qubit attached to another edge well. Unlike conventional 
approaches involving correlated voltage gate switches to operate on electron charge state, 
we make use of A - type resonant optical excitation scheme proposed earlier to drive 
the single-qubit evolution in various double- well quantum systems [7j - [21]. At most 
two resonant terahertz pulses are required in this case, and very accurate population 
transfer can be attained in times of order of several tens of picoseconds and for the 
pulses of much more moderate strengths than in the strategies mentioned above. Here 
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we study a way how to choose the pulse and structure parameters in order to optimize 
the efficiency of addressed probe electron transfer. We analyze the transfer probabilities 
at different choices of the pulse strengths and frequencies as well as for different numbers 
of quantum wells forming the structure. Although the numerical results are obtained 
for an effective one-dimensional structure, the algorithm of the structure and pulse 
engineering developed here can be adopted for realistic solid-state nanostructures. 

The paper is organized as follows. In Sec. II the parameters of individual quantum 
wells are selected in a way to form the reliable single-electron transport channel in the 
nanostructure composed of those wells. The energy spectrum and optical properties 
of multiwell nanostructure are considered in Section III. Section IV contains numerical 
data of the dynamical simulations. We conclude our study by Section V. 



2. Single-well parameter choice 

We begin with the description of the model of nanostructure employed in this work. The 
nanostructure is composed of quantum wells (QWs) placed in series and separated by 
finite barriers. Two edge QWs with indices i = 1 and i = N are deeper than the internal 
QWs. The energy spectrum of an electron confined in the structure is presented by two 
ground states (each of them being localized in corresponding edge QW) and N excited 
eigenstates delocalized over the structure. The latter result from the hybridization of the 
states of individual QWs due to the tunneling process. The locahzed states are supposed 
to be isolated from each other, i. e., the electron tunneling between them is negligibly 
small. If the electron occupies one of those states, it can stay there for extremely long 
time. The excited states of the nanostructure are approximately expressed by linear 
combinations of the states of isolated QWs [Ij. Provided that the nanostructure is 
symmetric, the wave functions of the hybridized states have even or odd parity relative 
to the center of nanostructure. The degree of hybridization depends on the heights 
and the widths of the potential barriers separating neighboring QWs as well as on the 
differences in the corresponding eigenenergies of isolated QWs. 

At the initial step, we choose the parameters of individual QWs. Several 
circumstances should be taken into account here. First, the ground states of the 
edge QWs are supposed to be strongly localized whereas the excited states of those 
wells must be close to the top of the potential barrier in order to form the reliable 
transport channel (through their hybridization with the states of neighboring internal 
QWs). Second, the matrix elements of optical dipole transitions (ODTs) between the 
localized and excited states of each edge QW have to be large enough to allow rapid 
optical excitation/deexcitation of an electron. It is known that the value of the ODT 
matrix element \do,k\ = |(0|ea;|fc)| (e is the electron charge, x is the coordinate) between 
the ground state |0) and an arbitrary excited state \k) of a single QW goes down quickly 
with the quantum number k of the excited state. For a square QW with infinite walls it 
depends on the excited state number k as |(io,fc| ~ (A; + 1)^ k'^ {k + 2)^ , where k is odd 
because of the selection rule. It means that each edge QW should be fabricated as a deep 
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narrow well containing only two bound electron states with /c = and k = 1. Further, 
since the tunnel coupling between edge QWs is mediated by an array of identical internal 
QWs, to implement electronic transfer over longer distances, one has to utilize wider 
QWs and to increase the QW number A^. From the other hand, it is desirable that 
each internal QW contain only one bound state. As a result, in order to meet above 
conditions the internal QWs should be designed as wide shallow wells. 

Instead of commonly used square potential shape, we consider here a less-studied 
power-exponential (PE) form for QW potentials |22j : 




where Uq is the depth of potential well, is the radius of potential, which determines 
the QW size, xq is the coordinate of the center of potential, and parameter p defines 
the smoothness of QW boundaries. This kind of electron confinement is found to be 
a good approximation for the electrostatically-gated quantum dots in two- and three- 
dimensional cases |22]. Wide class of potentials can be modelled within this framework; 
in particular, a square QW potential is obtained from Eq. (1) if one sets p >> 1. Note, 
that PE potentials are smooth functions of coordinate that considerably simplifies the 
numerical treatment of the eigenvalue problem. 

The one- dimensional Schrodinger equation for a single electron moving in the 
potential ?7(x). 



is solved with the help of the finite-difference second-order scheme [23]. This procedure 

is equivalent to the solution of a system of uniform algebraic equations expressed by five- 
diagonal sparse matrix and, consequently, to finding the eigenvalues and eigenvectors of 
that matrix. The accuracy of the method is very high, and the finite-difference schemes 
of higher orders are generally excessive for treating the one- dimensional problems. The 
results of calculations are given in effective atomic (or donor) units - i. e., 1 a.u. = 
Ry* = nfRy/meS^ for the energy and 1 a.u. = o*b = rneSaB/m* for the length, where 
Ry is the Rydberg energy, is the Bohr radius, mg is the free electron mass, m* is the 
effective electron mass, and e is the dielectric constant. For GaAs one has Ry* = 0.006 
eV and a^=10 nm. 

Here we demonstrate the solutions of Eq. (2) for a single QW with the PE potential 
of Eq. (1), where the potential radius and the parameter p are fixed and the potential 
depth Uq is varied. We used the computational interval of length L=10 with the number 
of grid points A'^2:= 10000 (at the boundaries of this interval the wave functions are set to 
zero). We take ao = 0.4 for the edge QW, ao = 1 for the internal QW, and p = 5, xo=0 
for both types of wells. The dependencies of low-lying part of the energy spectrum for 
the edge and internal QWs on Uq are presented in Figs. 1 (a) and 1 (b), respectively. 
We present the results for the values of Uq at which there are no more than two bound 
electron eigenstates with energies Eq < and ei < in the QW. Other states, whose 





(1) 



2m* dx'^ 



(x) + f/ (x) ^ (x) = (x) , 



(2) 
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energies are positive (hereafter - unbound or quasicontinuum states), belong to the 
spectrum of large square infinite potential well of the length L, modified by the QW in 
its central part. As the absolute value \Uq\ of the QW depth decreases, the levels Sq and 
El shift to the barrier top. They interact with unbound states of the same symmetry 
giving rise to the anticrossing pattern. The presence of the symmetric QW potential in 
the center of computation interval amounts also to the energy splitting of states of the 
large well. We observe apparent doublet structures for the states with positive energies, 
which transform into the states of the large square well at vanishing Uq. 

Figures 2 (a) and 2 (b) illustrate the dependencies of absolute values of ODT matrix 
elements \do^k\ on the potential depth Uq. The values of the ODT matrix elements are 
given in atomic units |e|a^. We see that the plot of |(io,i| for the edge QW [Fig. 2 (a)] has 
a local maximum at Uq ~ —24. At that point, other matrix elements of ODT are still 
small, and all excitations from the ground state of the edge QW to the quasicontinuum 
states remain inefficient. It is thus reasonable to take Uq ~ -24 for the edge QW. As 
the QW becomes shallower, one observes a sharp increase of the ODT matrix elements, 
which approach the values of the ODT matrix elements SqI of the large square well. 
The latter are much greater than the matrix elements of ODT between the bound 
states of deep QW; the values for both types of matrix elements relate each other as 
the corresponding effective potential lengths, namely, \d\^l/dQ^k\ ~ L/a^^ where L/oo 3> 
1. At very large potential depths the matrix element (io,i becomes close to the value 
(io°| = 16e/o/97r^ for an infinite square well, where Iq = 2ao is the well width, and for 
ao=0.4 one has d'^'l = 0.14. In contrast to the edge QW, all of the ODT matrix elements 
of internal QWs should be selected as small as possible in order to minimize unwanted 
excitations of electron from bound state to the quasicontinuum states. As it is seen in 
Fig. 2 (b), the plot of |(io,il demonstrates a pronounced increase for Uq > -2.5. Hence 
we choose -3 < C/q < -2.5 for the internal QW so that the first excited state is outside 
the QW but the values of the ODT matrix elements are still minor. 

Making our choice of the value of Uq for edge QWs [Fig. 1 (a)], we fix the electronic 
eigenenergies Eq and Ei of the edge QWs; then from the plot of the ground state energy 
Eq of internal QW vs Uq [Fig. 1 (b)] we define the value of Uo at which the energy 
£0 of internal QWs is approximately equal to the excited state energy Ei of the edge 
QWs. The results of calculations presented below are obtained for the following set of 
the parameters: oq =0.4, Uq = - 23.83 for the edge QWs and oq =1, Uq — - 2.5 for the 
internal QWs. From Figs. 1 (a) and 1 (b) one may see that the energy Si of the edge 
QWs sfightly differs from the energy Sq of the internal QWs. This is done intentionally 
in order to make our model of nanostructure more realistic since in practice it is very 
difficult to attain the exact adjustment of the QW levels. 

3. The nanostructure 

The N QWs stacked in the way described above constitute the body of a one-dimensional 
linear nanostructure. The energy spectrum of the nanostructure depends on the 
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parameters of individual QWs [the QW potential depths ?7o(0 ^^^^ the QW radii ao(^)) 
i = 1 - N] and on the widths of the barriers separating the wells. For simplicity, we 
suppose that all of the distances 6(2), i = 1 - {N — 1), between the centers of neighboring 
QWs (hereafter - interwell distances) are the same and equal to h. In this case, the 
potential of nanostructure is given by the expression 



where the origin of coordinate is placed into the center of nanostructure. Now we seek 
the solution of Eq. (1) with the potential having the form of Eq. (3) with the values of 
the edge and internal QW parameters selected above and p=5. The parameters for the 
edge (internal) QW are supplied with index (1): Uoii) = f/o(i), cto(0 = ^o(i) if i = l,N 



In what follows, we shall be interested in eigenenergies of not only the states bound 
in the QW structure, but also of states lying higher than the potential barrier. The 
accurate calculation of the energies of those quasicontinuum states requires a quite 
large computational length L and a great number A^^ of grid points. We have computed 
the eigenenergies Sk in the range Eq < Ek ^ \eo\ for different QW numbers N (up to 
= 20) at several sets of L and Nl- It was found that for all those nanostructures the 
choice L = 500 and Ni = 40000 enables one to satisfactory simulate the quasicontinuum 
states in the pointed energy interval. 

In this Section, we examine in detail the nanostructure composed of six QWs, see 
Fig. 3 (a). The plots of the excited state energies Ek {k=2 - 7) on the interwell distance 
h (actually, on the barrier width) are presented in Fig. 3 (b). As expected, the energies 
Ek depend on h exponentially. At the same time, they deviate from the tight-binding 
(TB) predictions [1] and bring about an asymmetrical distribution of Ek relative to the 
center of the energy subband. This is the peculiarity of a real system, where the energy 
of a hybridized state is to be a function of the closeness of this state to the barrier 
top. In fact, the overall correspondence between the spectrum of our structure and the 
spectrum of structure modelled within the TB approximation [IJ is mainly determined 
by the effective tunneling rates between neighboring QWs. Usually, it is expected that 
the weaker the tunneling the better this correspondence. Therefore, for deeper wells 
and wider barriers the results of TB model will give more accurate description of the 
nanostructure. Besides, the small detuning in the energies of individual QWs affects, to 
some extent, the mechanism of the state hybridization. For a quite large h {b >4.5), it 
results in a dissociation of the six-fold energy subband into the two degenerate energy 
levels {e2 and £3), presented by the even and odd superpositions of excited states of 
isolated edge QWs, and the four-fold subband of weakly- hybridized states with Ek^2,3 
pertaining to the internal part of the nanostructure. As follows from our calculations, 
the variation of the distance bo between the centers of the edge QWs and the neighboring 
internal QWs at some fixed distance bi between the centers of two neighboring internal 
QWs does not change noticeably the properties of the excited subband. The wave 
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functions of the bound states \k) (A;=0 - 7) are shown in Fig. 4. Note, that the 

wave functions '^o{x) and '^i{x) forming the ground state subspace of the Hilbert space 
of multiwell nanostructure are taken as the functions each localized in corresponding 
edge QW. This representation given by the isolated edge QW orbitals is equivalent to 
the representation given by their symmetric and antisymmetric combinations provided 
that the corresponding eigenenergies Eq and ei of the Hamiltonian are almost degenerate. 

Other important quantities, the ODT matrix elements, are plotted in Figs. 5 (a) 
and 5 (b). Due to the inversion symmetry of our structure, we may calculate the values 
of the matrix elements for only one excitation arm of the A scheme, e. g., involving the 
ground state |0), localized in the left QW («=1), and the excited states \k) {k=2 - 7). 
The ODT matrix elements for another excitation arm with the ground state |1) of the 
right QW {i=N), are related to them by the formula di^k = (— 1)'^'^o,A:- In the range of 
the interwell distances b around b =3, the values of matrix elements \do,k \ for pairs of the 
states (doublets) {|2), |7)}, {|3), |6)}, and {|4), |5)} remain close to the corresponding 
values calculated in the TB model (denoted by horizontal dotted lines). As in the case 
of the energies, both data types do not coincide. This discrepancy can be explained in 
the following way. The TB wave functions of hybridized states are constructed from the 
isolated QW wave functions multiplied by the weight coefficients that depend on both 
the state number k and the QW index i [1]. Consequently, the ODT matrix element 
between the ground state |0) and the hybridized state \k) is simply the single- well matrix 
element (io,i [Fig- 2 (a)] multiphed by the weight coefficient for the state \k) in the edge 
QW (^=1, A^). For two quantum states of a given doublet, those coefficients are the 
same in one of the edge QWs and differ by the sign in another one. Hence, the moduli of 
the ODT matrix elements |(io{i),fc| for doublet states turn out to be equal to each other 
and do not depend on b under the TB approximation. However, we observe in Fig. 4 
that the wave functions of doublet states have close but still different amplitudes in the 
edge QW regions. As a result, the matrix elements of ODT are different for the two 
doublet states. It means also that the wave functions of highly-delocalized states of the 
realistic nanostructure cannot be satisfactory approximated by the linear combinations 
of the wave functions of isolated QWs. 

If the QWs are substantially remote from each other (6 > 4.5), the curves for 
\do,2\ and |(io,3| demonstrate asymptotical convergence to a common non-zero value 
whereas the others decay slowly to zero. Such a behavior of the ODT matrix elements 
indicates on the collapse of tunneling between the edge QWs and the internal part of 
nanostructure. Since the states |2) and |3) evolve with b into even and odd superpositions 
of excited states of the isolated edge QWs, the linear combination {\do,2\ + Mo,3|)/"\/2 
tends to the value of |(io,i| for isolated edge QW. Note, that all of \do^k\, involving 
the quasicontinuum states \k) {k >7), are very small in comparison with the values 
of matrix elements for ODTs between the bound states. [Several matrix elements for 
the transitions from the ground state into the quasicontinuum states are shown as thin 
black curves at the bottom of Fig. 5 (a)]. The explanation of this fact is rather simple. 
In general, the one- dimensional wave function of a localized electron may be presented 
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in a form {x) = (x) j ^/h, where k is the quantum number of the electron state, 
Ik is the characteristic length of electron spreading in this state, and (x) is the 
dimensionless function which amplitude is of order of unity. The ODT matrix element 



between the states |0) and \k) can thus be written as (io,fe = Co,fce/Q / ^/loh, where 
IC'o.fcl <1- Substituting here Iq = Ik = ao, we see that corresponding matrix element of 
the ODT, involving two bound states |0) and \k), is of the order of eao- If, however, 
one sets Iq = oq and Ik = L, the estimation \do^k\ ~ gclq^Jo-o/L -C eao is obtained. Both 
results are consistent with the numerical data presented in Fig. 5 (a). 

Considering the next type of the matrix elements, we deal with the transitions 
between the states that constitute one subband. In Fig. 5 (b), we show the data 
illustrating the dependencies of such intrasubband ODT matrix elements on interwell 
distance b for the subbands {|0), |1)} and {|2) - |7)}. In contrast to the intersubband 
ODT matrix elements do^k {k = 2 -7) that vary smoothly with b, the intrasubband matrix 
elements are approximately linear functions of the interwell distance. We observe also a 
strong increase in the values of the intrasubband ODT matrix elements in comparison 
with the intersubband ones. Besides, the intrasubband ODT matrix elements between 
neighboring states are significantly larger than those between non-neighboring states, 
pictured as thin black curves in Fig. 5 (b). According to the TB model, all intrasubband 
matrix elements are expressed by linear combinations of the diagonal single-well matrix 
elements. The latter describe the shifts of the energy levels of isolated QWs in an 
electric field. Their difference has a clear analogy in the classical physics: it is just the 
energy acquired by an electron transported between the QW centers in the electric field 
of unit strength. In the case of the diagonal matrix element for the localized state |0), 
the "classical" formula |(io,o| = b{N — l)/2 perfectly reproduces the plot for |(io,o| vs b 
found in the essentially quantum treatment of the problem. With that, as b increases, 
the plots of matrix elements for excited subband substantially deviate from the linear 
form. For example, the matrix element (^2,3 initially separates from the others and then 
approaches (io,o at 6 > 4.5. Again, it is because of the reorganization in the structure of 
hybridized states due to the tunneling collapse. Thus, the intrasubband ODT matrix 
elements not only account of the effects associated with the charge displacement in the 
external field, but also reflect the peculiarities of electronic density distribution in given 
quantum states. 

The eigenenergies and matrix elements dm,n of ODT calculated above will be 
used in the next Section in numerical study of the dynamical problem. 

4. Dynamics 

We have mentioned in Introduction that there exist several ways of how to transport 
an electron from one QW to another in a semiconductor nanostructure. The aim of the 
present paper is to clarify some new aspects of resonant optical excitations in the effective 
A system formed in the spatially extended multiwell structure by a pair of ground states 
and a set of delocalized states. In the work [7], the symmetric three- level scheme was 
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employed to describe the resonant electron dynamics in the double-dot structure. It was 
shown that the square harmonic laser pulse of strength E (t) = Eq cos {ut) and duration 
T can produce complete population transfer between two spatially remote ground states 
|0) and |1) of the nanostructure if the pulse frequency u is tuned exactly to the resonance 
with the transitions {|0) , |1)} ^ |r) connecting the ground state subspace and some 
excited (transport) state |r) (other excited states are ignored). The accuracy of such a 
three-level approximation is very high provided that the following conditions imposed on 
the structure and pulse parameters are satisfied, (i) The transition frequency u is much 
larger than both the interlevel spacings |Ar ,,-|-i| = — e^-i-il and the coupling coefficient 
Ao,r = -Eo(io,r/2. It allows one to apply the so-called rotating wave approximation (RWA) 
under which an analytic solution of the dynamical problem is easily found, (ii) The 
coupling coefficient Ao,r is rather small as compared with |Ar,r±i|- It guarantees the 
selectivity of resonant excitations, i. e., only the transport state |r) is populated during 
the pulse action, (iii) The quantum dots (wells) are close enough to each other to 
neglect the diagonal coupling coefficients. It means that rapidly varying energy shifts 
do not modify appreciably the transition frequency co'o(i),r and thus do not violate the 
resonant conditions. The correction to the bare three- level scheme [7], accounting for a 
non-resonant excitation of the delocalized state nearest to the transport state, as well 
as the first-order correction to the RWA, were derived in Ref. [TO] . 

Here we investigate the problem of the resonant field-structure interaction via fully- 
numerical treatment. We solve the one-electron non- stationary Schrodinger equation 

zn^^\^{t)) = H{t)\^{t)) (4) 

with the Hamiltonian 

H (t) = Ho- cEoX cos (Lut) , (5) 

where Hq is the Hamiltonian of unperturbed nanostructure with the potential relief 
given by Eq. (3). 

The state vector of the system may be presented in the form 

I^W) = Ecn(t)e-^^"*/^|n), (6) 

n 

where index n runs over the states belonging to the low-energy part of the nanostructure 
spectrum. Inserting Eqs. (5) and (6) into Eq. (4) and then multiplying Eq. (4) from 
the left side by (m| (m = 0, 1, 2,..), we arrive at the set of linear differential equations 
for the probability amplitudes Cm {t). It is convenient to rewrite this set in the following 
way: 

.dc„ 



. 2 5]c„A™,„cos(cut)e-*"-"* (m = 0, 1, 2, ..), (7) 

where ujm,n = ^m,n/^ = i^n — £^m)/^ IS the transition frequency between the states 
|m) and \n), \m,n = £fieiddm,n/'^ is the corresponding coupling coefficient, Sfieid = 
ea*^EQ/Ry* is the field energy, and the time t is taken in units of h/Ry* (for GaAs 1 
a. u. =0.11 ps). 
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For definiteness, assume that at t=Q the electron is locahzed in the state |0) in 
the left edge QW - i. e., |\E' (0)) = |0), and hence the initial conditions for Eq. (7) 
read Cm (0) = 5mfl- Since our attention is concentrated on the possibility of complete 
population transfer between the ground states |0) and |1) of the multiwell nanostructure, 
here we restrict ourselves by analysis of the probability pi(T) = |ci(T)p to find the 
electron in the target state |1) after the pulse is switched off. The numerical solutions 
of Eq. (7) have been obtained with the help of the non-stiff odell3 Matlab solver with 
relative accuracy 6toi = 10~®. In all simulations we take the pulse frequency u to be 
equal to one of the nanostructure transition frequencies (-1^0(1),^ (2 < r < 7 for six- well 
structure). In this case, the single-electron dynamics is given by the resonant three- 
level Rabi oscillation picture. Typical plots of probabilities PoiT) and Pi(T) of finding an 
electron in the states |0) and |1), together with the total probability ptr (T) = J2 Pk (T) 

[pfc(T) = |cfe(T)p] of electron to be out of ground-state subspace, versus the pulse 
duration T, are presented in Fig. 6 for r = 5, = 6. 

The number of states, included into the simulations, was varied from 30 to 40. 
Apart from the bound states shown in Fig. 3 (a), we also took into account two groups 
of quasicontinuum states. The first group contains the states with energies < |£o| 
lying in nearest neighborhood of the barrier top. It is expected that those states may 
affect, to some extent, the resonant dynamics in the bound state subspace through non- 
resonant excitations. The second group incorporates the states from the narrow energy 
interval located around the energy Sr + uj. Since those states are close to the two-photon 
resonance with the driving pulse, they could participate the electron evolution despite 
of small values of the corresponding ODT matrix elements. However, as one observes in 
Fig. 6, the resonant dynamics of an electron in the multilevel system is an essentially 
three-level dynamics. It is very close to that considered in Ref. [7j and is described by 
the formulae 

Po = cos'^ (flRt) , pi = sin^ (flRt) , Pr = - sin^ (2fii?t) , (8) 

where Qr = |Ao,r|/ a/2 is the Rabi frequency. It means that only the states |0), and 
|r) play a decisive role in the electron dynamics, whereas the presence of other states is 
manifested in the weak modulations of ideal three- level plots. The participation of ladder 
resonant excitations involving quasicontinuum states as well as the processes beyond the 
RWA may be safely ignored due to the smallness of the parameters |Ar, A;|/|Ao,r| {k > 7) 
and |Ao,r|/^- The transfer time Tq is thus defined from the condition VL^Tn = 7r/2 + vrn 
(n=0, 1,2,..) as the shortest pulse duration T!„ at which the maximal population of 
target state |1) is achieved - i. e., Tq = 7r/2f2^. Note, that the definition of Rabi 
frequency ^Ir is taken according to Refs. [7] and [10]; it differs by 1/V2 from that used 
in Ref. [llj for the description of resonant hydrogen molecular ion dynamics [cf. |llj . 
Eq. (40)], since in that work the ground states of a molecular ion, as well as its excited 
states, were presented by even and odd linear combinations of the isolated atom states. 

The most important question one has to answer is what a state from the excited 
subband should be used as the transport state |r) in order to optimize the transfer 
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probability max(pi) and to reduce the transfer time Tq. According to the selectivity 
and high speed requirements, the values of |(io,r| and |A.r^r±i| should be taken as large as 
possible. As it is clearly seen from both numerical data and analytical TB estimations, 
the states from the center of the excited energy subband (in the six-well structure the 
states |4) and |5)) are characterized by the largest interlevel spacings and, at the same 
time, by the highest ODT matrix elements |(io,r| than the states at the subband edges. 
We have performed the calculations of maximal transfer probabilities and transfer times 
at several fixed field energies e field (actually, the field strength) and for the different field 
frequencies uo = cuq,*: {k=2 - 7). The results of calculations confirm our expectations, 
demonstrating that To,fc=4 < To^k=5 < ^0,^=3 < To,fc=6 in the whole range of e field- The 
pulse durations TQ k=2 and TQ^k=7 corresponding to the electron transfer via edge excited 
subband states are significantly longer than others. Thus we come to the conclusion that 
the state closest to the center of excited subband is the best candidate for the transport 
state. If the QW number N is even, it will be one of states {\N/2 + 1), \N/2 + 2)} of 
the central doublet, whereas for odd this is the unpaired state |1 + (A^ + l)/2). 

The plots of the maximum values of pi, see Fig. 7, demonstrate quasiperiodic 
behavior of the transfer probabilities on the field energy for each state of the pair {|4), 
1 5)} for A^=6. We observe that the state |4) displays better transport characteristics 
in comparison with the state |5). The maximum of the transfer probability max(]9i) ^ 
0.999 for = Ci;o,4 is attained at quite large field energy e field = 0.083, thus allowing for 
rapid electron transfer (about 45 ps for the GaAs parameters). The transfer times Tq for 
both = ci;o,4 and uj = ujq^^ cases are presented in Fig. 8. It is worth mentioning that 
the simple three-level expression for Tq (shown by the dashed curves) fits the numerical 
graphs very well. 

In order to explain the dependencies of max(pi) on the field energy, we make 
a natural assumption that the oscillations in corresponding plots are caused by the 
presence of other excited states close to the transport state |r). In the work [10], we 
calculated the correction to ideal resonant three-level dynamics [Eq. (8)] taking into 
account the non-resonant occupation of the neighboring state \k). The corrected four- 
level formula for the population pi of target state |1) at T = Tq reads in current notations 

(To) ^ 1 - /,,, = (2Ao,fc/Afc,,)'sin2 (vrA,, ,/2v^Ao,.) . (9) 

Note that this formula was derived via accurate solution of the four-level dynamical 
problem and contains more profound physics than those obtained by the simple adiabatic 
exclusion of the non-resonant levels. 

The form of the set of Eqs. (7) assures that the non-resonant excitations of the 
delocalized nanostructure states \ k) {k^r) occur independently up to the terms of higher 
orders in the parameters Xo^k/^k,r, and the direct generalization of Eq. (9) to the 
multilevel A system results in the formula 

max(pi)=pi(To)^l-^/fc (10) 
(hereafter, the sum is over the bound states only). The dependencies of max(pi). 
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calculated through the Eq. (10), are visualized in Fig. 7 as the dashed curves. We see 
that for u — a;o,4 the formula (10) is in a good agreement with the numerical data, while 
for uj = cjo.s it gives overestimated values for local maxima of the transfer probability. 
What is the reason for such deviations? To our opinion, the main error originates from 
the neglecting in fk = |cfcp [Eqs. (9) and (10)] the interference terms describing the 
correlation effects between the state \k) and other excited states. Despite of the small 
values of these terms in comparison with non-correlated terms entering the formula 
(10), the correlated excitations may amount to some reduction of transfer probability. 
Besides, in the derivation of Eq. (9), the diagonal (intrasubband) coupling coefficients 
were ignored. These intrasubband terms of Eq. (7) oscillate at the field frequency cu 
relative to the RWA frame and are usually discarded. However, since their values grow 
linearly with the interwell distance, they will affect the electron dynamics in extended 
structures, as \m,n become close to u. 

We believe that the difference in the influence of one excited state on another 
excited state is explained by the difference in their ODT matrix elements |rfo,fc| • For 
example, the value of \do^^\ is larger than the value of |(io,5|; as a result, the population 
of the state |5) is subjected to stronger modulations caused by the state |4), while the 
influence of the state |5) on the the population of the state |4) is less pronounced. It 
means that the state |4) possesses better transport properties than the state |5) and thus 
it should be selected as the transport state in the six-well nanostructurc. The origin of 
oscillations of the transfer time Tq may be interpreted in a similar way. It should be 
noted that the first-order oscillating correction to Tq docs not improve the accuracy of 
Eq. (10); this observation agrees with the arguments given above that several processes 
were neglected in derivation of Eq. (10). Nevertheless, the analytical curves yield rather 
accurate values of the local maxima positions in both cases. Therefore, Eq. (10) may 
be used at initial step for evaluation of the laser field strengths for which the probability 
of electron transfer is maximal. On the other hand, this formula enables us to study 
the dependence of max(pi) on the energy splittings that, in their turn, are the 
functions of the interwell distance b and the QW number N. 

The important feature in the behavior of max(pi), expressed by Eq. (10) and 
visualized in Fig. 9 for = 6, is its band-like structure, where high transfer probability 
regions alternate with low transfer probability regions. We observe that max(pi) exhibits 
oscillations as b varies. Numerical data, plotted in Fig. 7 (a) for 6 = 3, lie at the vertical 
dashed fine in Fig. 9. We have found that the value of max(pi)=0.999 at e field— 0-0S3 
(the point ^4) is likely to be the optimal one, since the closest maxima of max(pi) are 
lower. As predicted by the formula (10), if ejieid approaches 0.14, the dashed line crosses 
the region where the transfer probability substantially increases. At the point B (6=3, 
5yjeZd=0.14) the analytical value of max(pi) = 0.991 correlates well with the numerical 
value of max(pi) = 0.988. It should be noted that, despite of the appreciable decrease in 
the transfer probability in comparison with the point A, the value of transfer time Tq = 
250 for e field — 0.14 is about half of Tq = 420 for e field — 0.083. We have also performed 
the calculations of max(pi) for the set of field and structure parameters chosen from the 
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high probabihty band containing the points C and D (see Fig. 9). Again, the numerical 
data are close to the analytical results, and max(pi) grows monotonically from 0.94 to 
0.98, as we move from C to D along the dashed broken line, whereas the transfer time 
decreases from 419 to 252. Going beyond the point D, we arrive at the crossing region 
between this band and the vertical dashed line marking 6 = 3, that covers the field 
energy interval around e field = 0.2 (not shown). For e field = 0.2 we obtain max(pi) = 
0.986 that is smaller than the value computed at the point B, but the corresponding 
transfer time Tq = 175 appears to be very short. It means that if we are interested in 
fast electron transfer with moderate reliability, it is recommended to take e field = 0.2 
instead of e field = 0.083. At large values of the parameters e field and b, the probability 
of successful electron transfer becomes quite low (dark area in a right-up part of Fig. 
9). This result is consistent with that obtained earlier (see Refs. |T0] and [H]) and is 
explained in terms of the tunneling collapse. 

At the final step, we have performed numerical simulations on the electron dynamics 
in the nanostructure composed of N QWs, where the QW number varies from = 5 
to A^ = 20, and 6 = 3. If A^ is not large (5 < A^ < 12), the dependencies of max(pi) 
on the field energy e field resemble those shown in Fig. 7 and demonstrate quasiregular 
oscillations with the amplitude increasing smoothly with A^. For A^ > 17, the plots of 
max(pi) versus e field exhibit a "sawtooth" profile with sharp deeps and peaks. Again, 
we make use of Eq. (10) as the approximation for max(pi), when the central state of the 
excited subband serves as the transport state. We find that the shapes of numerical and 
analytical curves of max(pi) deviate substantially from each other, if A^ approaches 20, 
whereas the correspondence between the values of peak positions remains surprisingly 
good. A general tendency is the shrinking of the region, where the values of max(pi) 
exceed 0.99. The last pronounced peak of max(pi) moves from e field = 0.2 for A^ = 6 
to e field = 0.08 for A^ = 20. We reveal, however, that the transfer times Tq taken at 
the peaks of max(pi) are very close to each other for different A^. In Fig. 10 we present 
the dependencies of Tq on e field for A^ =15, 16, 17, and 18. The points of interest, 
where 0.99 < max(pi) < 0.999, are marked by open (filled) squares for odd (even) A^. 
It is easy to see that each group of peaks, pertaining to different A^, corresponds to a 
narrow range of Tq values. This observation enables us to make an important conclusion 
that the electron transfer time in multiwell nanostructure is not strictly defined by the 
number of QW in the structure, as it might be expected, but displays rather complex 
behavior conditioned by both the structure and field parameters. As the final remark, 
we note that the optimal points of the transfer probability lie in the middles of graph 
steps, where the analytic expression Tq = 7r/|Ao,r|"\/2 yields very accurate values of the 
electron transfer time. Thus it is expected that this simple formula can be used to 
estimate Tq in the multiwell nanostructures regardless the QW number A^. 
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5. Conclusions 

In this paper, we have investigated theoretically the important and very actual 
problem of coherent single-electron transfer between two remote quantum wells in the 
semiconductor linear nanostructure. Apart from wide interest to the fundamental 
question addressing the possibility of reliable control over individual charge carrier 
dynamics in the solid-state low-dimensional structures, this study has a clear practical 
application - i. e., the organization of indirect connections in the ordered arrays of 
quantum bits. 

For this purpose we make use of an auxiliary nanostructure containing a single 
electron in the quantized part of its conduction band. The resonant monochromatic 
laser pulse, applied to the structure, drives the electron evolution in such a way that, 
being initially localized in one quantum well, it comes onto another well at the end of 
the pulse. This process is influenced by many factors. In order to optimize the protocol 
- i. e., to attain high probability of electron transfer in short time, one should to choose 
carefully both structure and external field parameters. As we have found, the maximum 
electron transfer probability is a quasiperiodic function of the pulse strength and the 
interlevel spacings in the excited subband. The choice of one of excited states close to 
the subband center as the transport state is substantiated by an observation that this 
state possesses large ODT matrix element and is far enough from the others. It provides 
high transfer characteristics such as speed and selectivity. 

Finally, it may be expected that the driving schemes relied upon the stimulated 
Raman adiabatic passage (STIRAP) |T71 |2ll [25] together with optimal control 
techniques [20] will enhance the efficiency of optically-induced single-electron transfer 
in the multiwell nanostructures. 
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Figure 1. (color online). The dependencies of electronic eigenenergies £k of low- lying 
cigcnstates of a single PE quantum well on the potential depth Uq for a) the edge 
quantum well and b) the internal quantum well. The well parameters used in further 
numerical calculations are shown in the insets. 



Figure 2. (color online). The values |do,fe| of matrix elements of optical dipole 
transitions, connecting the ground state |0) and the low- lying excited states |A;) for 
a) the edge quantum well and b) the internal quantum well. 



Figure 3. (color online). The six- well one-dimensional nanostnicture placed in the 
central part of the large quantum well {L = 100 a.u.). a) Potential profile U{x) of the 
structure in the case of 6o = 6i = 3 a.u. and the energy levels Sk from ground {k=Q, 
1) and first-excited (fc = 2 - 7) subbands. b) The dependencies of the excited state 
energies Sk (fc = 2 - 7) on the interwell distance b — — bi. Vertical dotted line marks 
the value 6 = 3 a.u. at which the most of numerical simulations are performed. Other 
parameters are given in the plot. 



Figure 4. (color online). The plots of the electron wave functions for low- lying states 
of the six-well nanostructure (in arbitrary units) vs the coordinate x at b~3 a.u. a) 
Ground-state wave functions \l/o (x) and {x); b) Excited-state wave functions \E'fe (x), 
fc = 2 - 7. 



Figure 5. (color online). The absolute values of ODT matrix elements of the six- well 
nanostructure as the functions of interwell distance b. a) Intersubband matrix elements 
\do,k\ ■ The horizontal dotted lines designate the values of |do,fe| calculated in the TB 
approximation, b) Intrasubband matrix elements Irfm.nl- 



Figure 6. (color online). The populations and pi of the ground states |0) and |1), 
and the total population ptr of excited states states \k), k >2, vs the pulse duration 
T. The structure and pulse parameters as well as the value of the transfer probability 
max(pi) are shown on the plot. 
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Figure 7. (color online). The maximum transfer probability max(pi) of electron in 

the six- well nanostructurc plotted as a function of the field energy e field for two choices 
of the transport state: a) r=4 and b) r=5. The numerical solutions are presented by 
the solid curves (the calculated points are denoted by circles); the approximations of 
Eq. (10) are shown by the dashed curves. 



Figure 8. (color online). The transfer time Tq between the ground states of the 
six- well structure vs the field energy e fi^id for two choices of the transport state, r=4 
and r=5. The plots visualizing the approximation Tq = 7r/|Ao, r|\/2 are presented by 
dashed curves. 



Figure 9. (color online). The transfer probability max(pi), Eq. (10), as a function of 
both the interwell distance h and the field energy e field- The dashed lines denote the 
paths along which the numerical simulations were performed. 



Figure 10. (color online). The transfer times To vs the field energy e field for the 
QW numbers N =15, 17 (blue dashed curves, open squares) and =16, 18 (red solid 
curves, filled squares). The analytic data are shown by thin dotted curves. 



